EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



CERN-EP/2006-039 
Dec. 8 2006 

Measurements of Charged Kaon Semileptonic Decay 

Branching Fractions K —> and K — > and 

O ■ Their Ratio 

o 

J.R. Batley, C. Lazzeroni, D.J. Munday, M.W. Slater, S.A. Wotton 
Cavendish Laboratory, University of Cambridge, Cambridge, CB3 OHE, U.K. ^ 

R. Arcidiacono 2 ) , G. Bocquet, N. Cabibbo, A. Ceccucci, D. Cundy 3 ), V. Falaleev, 
M. Fidecaro, L. Gatignon, A. Gonidec, W. Kubischta, A. Norton 4 \ M. Patel, A. Peters 

CERN, CH-1211 Geneva 23, Switzerland 

S. Balev, P.L. Frabetti, E. Goudzovski, P. Hristov 5 ), V. Kekelidze 5 ** , V. Kozhuharov, 
L. Litov, D. Madigozhin, E. Marinova, N. Molokanova, I. Polenkevich, Yu. Potrebenikov, 

S. Stoynev 6 ), A. Zinchenko 
t^. ■ Joint Institute for Nuclear Research, Dubna, Russian Federation 

E. Monnier 7 ), E. Swallow, R. Winston 



(N 

(N 
> 

in 



o 



X 

D | The Enrico Fermi Institute, The University of Chicago, Chicago, IL 60126, U.S.A. 

gn! P. Rubin, A. Walker 

<-j ', Department of Physics and Astronomy, University of Edinburgh, 

JCMB King's Buildings, Mayheld Road, Edinburgh, EH9 3JZ, U.K. 

^ . W. Baldini, A. Cotta Ramusino, P. Dalpiaz, C. Damiani, M. Fiorini, A. Gianoli, 

M. Martini, F. Petrucci, M. Savrie, M. Scarpa, H. Wahl 
Dipartimento di Fisica dell'Universita e Sezione dellTNFN di Ferrara, 1-44100 Ferrara, Italy 

A. Bizzeti 8) , M. Calvetti, E. Celeghini, E. Iacopini, M. Lenti, F. Martelli 9) , 

G. Ruggiero 5) , M. Veltri 9 ) 
Dipartimento di Fisica dell'Universita e Sezione dell'INFN di Firenze, 1-50125 Firenze, Italy 

M. Behler, K. Eppard, K. Kleinknecht, P. Marouelli, L. Masetti, U. Moosbrugger, 
C. Morales Morales, B. Renk, M. Wache, R. Wanke, A. Winhart 
Institut fur Physik, Universitat Mainz, D-55099 Mainz, Germany 10 ^ 

D. Coward 11 ^, A. Dabrowski, T. Fonseca Martin 5 \ M. Shieh, M. Szleper, M. Velasco, 

M.D. Wood 12 ) 

Department of Physics and Astronomy, Northwestern University, 
Evanston, IL 60208-3112, U.S.A. 

G. Anzivino, P. Cenci, E. Imbergamo, M. Pepe, M.C. Petrucci, M. Piccini 5 ), M. Raggi, 

M. Valdata-Nappi 

Dipartimento di Fisica dell'Universita e Sezione dell'INFN di Perugia, 1-06100 Perugia, Italy 

C. Cerri, G. Collazuol, F. Costantini, L. DiLella, N. Doble, R. Fantechi, L. Fiorini 13 \ 
S. Giudici, G. Lamanna, I. Mannelli, A. Michetti, G. Pierazzini, M. Sozzi 



Dipartimento di Fisica, Scuola Normale Superiore e Sezione dell'INFN di Pisa, 

1-56100 Pisa, Italy 

B. Bloch-Devaux, C. Cheshkov 5) , J.B. Cheze, M. De Beer, J. Derre, G. Marel, 
E. Mazzucato, B. Peyaud, B. Vallage 
DSM/DAPNIA - CEA Saclay, F-91191 Gif-sur-Yvette, France 

M. Holder, A. Maier 5 ), M. Ziolkowski 
Fachbereich Physik, Universitat Siegen, D-57068 Siegen, Germany 14 ^ 

S. Bifani, C. Biino, N. Cartiglia, M. Clemencic 5 -* , S. Goy Lopez, F. Marchetto 
Dipartimento di Fisica Sperimentale dell'Universita e Sezione dell'INFN di Torino, 

1-10125 Torino, Italy 

H. Dibon, M. Jeitler, M. Markytan, I. Mikulec, G. Neuhofer, L. Widhalm 
Osterreichische Akademie der Wissenschaften, Institut fur Hochenergiephysik, 

A-10560 Wien, Austria 15 ^ 



Published in European Physical Journal C 50 2 (2007), contains erratum submitted to 

European Physical Journal C 



x ) Funded by the U.K. Particle Physics and Astronomy Research Council 

2 ) Present address: Dipartimento di Fisica Sperimentale deH'Universit Sezione dell'INFN di Torino, 
1-10125 Torino, Italy 

3 ) Present address: Istituto di Cosmogeofisica del CNR di Torino, 1-10133 Torino, Italy 

4 ) Present address: Dipartimento di Fisica dell'Universita e Sezione Ferrara, 1-44100 Ferrara, Italy 

5 ) Present address: CERN, CH-1211 Geneva 23, Switzerland 

6 ) Present address: Department of Physics and Astronomy, Northwestern University, Evanston, IL 60208- 
3112, U.S.A. 

7 ) Also at Centre de Physique des Particules de Marseille, IN2P3-CNRS, Universite de la Mediterranee, 
Marseille, France 

8 ) Also Dipartimento di Fisica, Universita di Modena, 1-41100 Modena, Italy 

9 ) Istituto di Fisica, Universita di Urbino, 1-61029 Urbino, Italy 

10 ) Funded by the German Federal Minister for Education and research under contract 05HK1UM1/1 
n ) Permanent address: SLAC, Stanford University, Menlo Park, CA 94025, U.S.A. 

12 ) Present address: UCLA, Los Angeles, CA 90024, U.S.A. 

13) p resen t address: Institut de Fisica d RAltes Energies, Universitat A Barcelona, E-08193 Bellaterra 
(Barcelona), Spain 

14 ) Funded by the German Federal Minister for Research and Technology (BMBF) under contract 056SI74 

15 ) Funded by the Austrian Ministry for Traffic and Research under the contract GZ 616.360/2-IV GZ 
616. 363/2- VIII, and by the Fonds fur Wissenschaft und Forschung FWF Nr. P08929-PHY 



Abstract 



Measured ratios of decay rates for 'R-Ke3lK2-Ki T^K\&Iki-k an d T^K^/KeS are P re_ 
sented. These measurements are based on K decays collected in a dedicated run 
in 2003 by the NA48/2 experiment at CERN. The results obtained are TZKe3/K2n = 
0.2470 ± 0.0009(stot) ± 0.0004(syst) and K K ^ 3/K2n = 0.1637 ± 0.0006(siat) ± 
0.0003(sysi). Using the PDG average for the — * 7r ± 7r° normalisation mode, 
both values are found to be larger than the current values given by the Particle 
Data Book and lead to a larger magnitude of the \V US \ CKM element than previ- 
ously accepted. When combined with the latest Particle Data Book value of \V u d\, 
the result is in agreement with unitarity of the CKM matrix. In addition, a new 
measured value oiTZ Kfl3 / Ke3 = 0.663 ±0.003(siai) ± O.OOl(sysi) is compared to the 
semi-empirical predictions based on the latest form factor measurements. 



1 Introduction 

New measurements of the charged kaon semileptonic decays, K — > 7r°e ± i> (K e3 ) 
and K — *> (K^), are presented. These measurements are based on K ± decays 

collected in 2003 by the NA48/2 Collaboration at the CERN SPS. 

The measured ratios are: 

_ r(K± -> TT°e ± u) _ r(K ± -> 77°/^) 

n K e3/K2n = _^ > T^K^/K2n = ~ I 1 ) 

and 

_ r(K ± -> ttV^) 

K W e3 = m± ~ ■ (2) 

In both the numerator and denominator, the decays contain a charged track and at least 
two photons originating from a 7r° decay, thus leading to a partial cancellation in the 
acceptance and reconstruction uncertainties. Contributions from internal bremsstrahlung 
are included for all three decay modes. The general reconstruction methods are the same 
for all three measurements but the event selection varies because different particle iden- 
tification criteria are applied. 

The main interest in measuring these quantities is to extract: (1) the individual 
semileptonic decay widths needed to determine the V us element in the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix; and (2) the ratio T{K^)/T{K e3 ) which is a func- 
tion of the slope parameters of the form factors. On the assumption of fi — e universality 
this ratio provides a consistency check between measurements of the form factors and of 
the partial decay widths. 

This paper is organized as follows. Section [2] outlines the phenomenological descrip- 
tion of semileptonic decays. The experimental setup is described in Section [3j The event 
selection and all corrections applied to the data are presented in Section HI Section 
includes a detailed description of the acceptance corrections and radiative effects. Finally, 
the results are presented in Section [HI and their comparison with theory and their impact 
on V us are given in Section 

2 Decay rates for semileptonic decays 

The decay rate for charged semileptonic decays can be written as follows [I]: 

W = ^<Sew\v us \ 2 \u(q)\ 2 i £ k (i + 4) , (3) 



where I refers to either e or /i, Gf is the Fermi constant, rriK is the kaon mass, Sew is the 
short-distance radiative correction, (1 + 5k) — (l + ^st/(2) ~^^em) 2 * s the model-dependent 
long-distance correction with contributions due to isospin breaking in strong (SU(2)) and 
electromagnetic (EM) interactions, /+(0) is the form factor at zero momentum transfer 
(t = 0) for the iv system. The remaining term, 1^, is the result of the phase space 
integration after factoring out /+(0), and is defined as pQ: 
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x {T(t,mLm2)(2t + m £ 2 )|/ + (t)// + (0)| 2 + 3mKm^-m2)Vo(t)// + (0)| 2 }, (4) 

where T(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz, and m n and nig are the masses of the 
7T° and the charged lepton, respectively. The two form factors correspond to the angular 
momentum configuration of the K — tt system, with f+(t) representing the vector form 
factor, while fo{t) is the scalar form factor. The t dependence of the form factors can be 
described using a quadratic (linear) approximation for the vector (scalar) term: 



f + (t) = f + (o) (l + \' + -^- + \kS~) and /o(t) = /+(0) ( 1 + Ao ^) 



(5) 



where /+(0) is obtained from theory, and A^, A'| and Ao are measured [2]. The second 
term in l l K shows that decays are more sensitive to the scalar form factor than K e ^. 

Other ways to describe the momentum dependence of the form factors are also 
considered and discussed in Section [5j 

Assuming /z — e universality, the Y^K^) /T(K e3 ) ratio is predicted to be [3]: 

0.645 + 2.087A+ + 1.464A + 3.375Aj + 2.573A 2 
TZ K , 3/Ke3 = T(K, 3 )/T(K e3 ) 1 + 3 .457A + + 4.783A 2 , ' (6) 

This semi-empirical formula assumes a linear approximation for the form factors, /+ o(t) = 
/ + (0)(l + A +i0 ^-). 

If the assumption of /i — e universality is removed, Eq. implies that: 

n K , 3/Ke3 = [g»mo)/get + m 2 x + ^))/(^(l + (7) 

where g is the weak coupling constant for the lepton current. The (1 + 5 e K ) /(l + 5^) ratio 
for charged kaons is very close to unity because the S^s are dominated by the few percent 
SU(2) correction that is in common between the electron and the muon channel. The EM 
corrections are at the per mille level and compatible with zero within errors [H El E] • This 
is not necessarily true for neutral kaons because there is no SU(2) correction and also the 
electromagnetic correction for the muons is larger. 

3 Experimental setup 
3.1 Beam 

The experiment uses simultaneous K + and K~ beams produced by 400 GeV pro- 
tons impinging on a Be target. The beam has particles of opposite charge with a central 
momentum of 60GeV/c and a momentum band of ±3.8% produced at zero angle. Both 



beams are selected by a system of dipole magnets forming an achromat, along with fo- 
cusing quadrupoles, muon sweepers and collimators. The spill length is about 4.5 s out 
of a 16.8 s cycle time, and the proton intensity is fairly constant during the spill with a 
mean of 5 x 10 10 protons per spill. The positive (negative) kaon flux at the entrance of 
the decay volume is 3.2 x 10 6 (1.8 x 10 6 ) particles per spill. 

For the measurements presented here, the proton beam intensity is reduced from its 
nominal value so that the data-acquisition system can handle the rate of the minimum bias 
trigger. The K + / K~ ~ 1.78 flux ratio is given by their production rate at the Be-target. 

3.2 Detector 

The final beam collimator is immediately followed by a 114 m long cylindrical tank, 
which is evacuated in order to minimize the interactions of the decay products. The tank 
is terminated by a 0.3% radiation length (Xq) thick Kevlar window, except in a region 
close to the beam which continues in a vacuum pipe through the centre of the downstream 
detectors. The detector is located downstream of this Kevlar window. 

The tracking is performed with a spectrometer housed in a helium gas volume. It 
consists of two drift chambers before and two after a dipole magnet with a horizontal 
transverse momentum kick of 120MeV/c. Each chamber has four views, each of which 
has two sense wire planes. The resulting space points are typically reconstructed with a 
resolution of 150 /im in each projection. The momentum resolution of the spectrometer is 
(Tp/p = 1.02% © 0.044% • p, where p is in GeV/c. The track time resolution is 1.4ns. 

The detection and measurement of electromagnetic showers are performed with a 
27 Ao-deep liquid krypton calorimeter (LKr) that has an energy resolution[7j of a(E)/E = 
3.2%/v / £©9%/£©0.42%, where E is in GeV. The calorimeter is subdivided into 13248 
cells of transverse dimension 2 cm x 2 cm. 

The reconstructed kaon mass in — ► tt^tt (K2tt) decays has a typical resolution 
of 3.25±0.05 MeV/c 2 , without imposing the 7r° mass constraint. 

A scintillator hodoscope is located between the spectrometer and the LKr for trig- 
gering purposes. It consists of two planes, segmented in horizontal and vertical strips 
respectively, with each plane arranged in four quadrants. The time resolution for the 
hodoscope system is 200 ps. 

Downstream of the LKr calorimeter is an iron-scintillator sandwich hadron calorime- 
ter (HAC), followed by muon counters (MUC) which consist of three planes of plastic 
scintillators shielded by the HAC and a 80 cm thick iron wall. The first two planes are 
made of 25 cm wide horizontal and 25 cm wide vertical scintillator strips, with a length of 
2.7m. The third plane consists of horizontal strips of width 44.6 cm and is mainly used 
to measure the efficiency of the counters in the first two planes. The central strip in each 
plane is split with a gap of 21 cm to accommodate the beam pipe. The fiducial volume of 
the experiment is principally determined by the LKr calorimeter acceptance. 

3.3 Trigger and readout 

The data coming from the detector is sampled every 25 ns, and eight samples are 
recorded in time windows of 200 ns. 

The trigger selection is defined by at least one hit in each of the horizontal and 
vertical planes of the hodoscope, within the same quadrant. The typical rate of this 
trigger at the reduced intensity is 1.3 x 10 5 per burst, and is downscaled by a factor of 
four in order to match the allowed limit of the data-acquisition system of 50000 events 
per burst. The efficiency of this trigger for events containing one track is calculated from 



events collected by requiring at least three hit wires in at least three views of the upstream 
chamber recorded within 200 ns. As shown in Table [TJ the trigger efficiency is found to be 
high and independent of the decay mode and the type of kaon beamline. 

4 Event selection and reconstruction 

The signals for the K e3 , K^ 3 and K 2lT decay modes require one track in the drift 
chambers and at least two clusters (photons) in the LKr that are consistent with a tt° 
decay. It is an important feature of this analysis that their basic event reconstruction and 
selection criteria are the same. These decay modes are separated from each other on the 
basis of particle identification and additional kinematic requirements, which also serves 
to reduce the background to a negligible level. 

4.1 Common selection criteria for K e3 , K^ 3 and K 27r 

All events are required to have at least one track that satisfies the following condi- 
tions: 

— the track must be at least 15 cm away from the center of each drift chambers and at 
least 15 cm from the centerer of the LKr, 11cm from its outside borders, and 2 cm 
from any inefficient cell; 

— to be within the allowed 37 ns readout window defined by the trigger hodoscopes; 

— to pass all track quality requirements. 

The decay vertex, V Xt y jZ , is reconstructed from the charged track and the beam axis. After 
correcting for residual magnetic fields inside the decay pipe, it is required that: 

— the longitudinal position, V z , must be 1316 cm downstream of the final beam colli- 
mator, and 2584 cm upstream from the end of the decay pipe, for a total of 75 m of 
decay volume; 

— the transverse position, V x>y , must be within three sigma (around 2.5 cm) in x and 
y from the beam axis. The actual cut depends on the longitudinal vertex position 
due to the varying uncertainty of the extrapolation to the vertex position. 

All cluster candidates in the 7r° reconstruction need to be clearly identified as photons, 
that is: 

— the energy of the cluster is greater than 3 GeV and less than 65 GeV; 

— time difference between clusters is smaller than 5 ns; 

— the minimum distance between clusters is 10 cm in order to minimize the effect of 
energy sharing on cluster reconstruction; 

— the minimum distance between clusters and extrapolated tracks is greater than 
10 (35) cm, so that no charged lepton (hadron) track is associated to the clusters. 

To make an inclusive measurement of the radiation, a minimum of two reconstructed 
photons is required. The threshold to define a reconstructed photon, 3 GeV, leads to a 
high fraction of the events with only two photons reconstructed. The fraction of events 
with at least one extra cluster in the calorimeter is 2.3% in the case of K e3 , 0.3% in 
and 10% in K2 n - Fig. [1] shows the invariant mass distribution for these two photon events 
reconstructed using the V z component of the charged vertex. The non-Gaussian tails are 
of the same order in all three channels. No requirement on the 77 invariant mass is made. 

For the small fraction of events with extra clusters, the best 7r° candidate is selected 
by comparing the V z component of the charged vertex and the longitudinal vertex, z n , 
of the 7T°. This neutral vertex is reconstructed by assuming the 7r° mass and using the 
kinematical information from the two photon candidates in the 7r° — > 77 decay. z v is 



calculated from the distance, d w , of the it vertex to the LKr calorimeter calculated from: 

d n = ^EiEj[(xi - xj) 2 + - yj) 2 ]/m n o , (8) 

where Xi and yi are the coordinates of the impact point of a the i th 7 in the LKr, and 
Ei is the energy of the corresponding 7. The longitudinal vertex position, z n , is given 
by z-n = z (LKr position) — d^. Although there is no explicit requirement on V z — z n , the 
combination with the smallest V z — is taken. 

4.2 Selection criteria to distinguish between K e3 , K^ 3 and K 2l , 

The invariant mass for the K 2lx candidates is required to be within three sigma of 
the reconstructed K ± mass, that is: 

0.4772 GeV/c 2 < m^o < 0.5102 GeV/c 2 , 

while for candidates: 

m^o < 0.4772 GeV/c 2 or m^o > 0.5102 GeV/c 2 , 

assuming m v ± instead of mi. The Monte Carlo simulated reconstructed invariant mass for 
the track and the ir°, without applying particle identification and under the assumption 
that the track is a 7r ± , is shown in Fig. [2] for the relevant decay channels. 

The allowed range in the transverse momentum of the track, Pp, with respect to 
the beam axis is required to be: 

K e3 : P T < 0.200 GeV/c, : P T < 0.200 GeV/c, K 2n : P T < 0.215 GeV/c. 

The missing mass squared, m 2 , is reconstructed under the assumption that the kaon 
energy is 60 GeV. The allowed m 2 range for each channel is: 

K e3 : - 0.012 GeV 2 /c 4 < m 2 < 0.012 GeV 2 /c 4 , 
: - 0.010 GeV 2 /c 4 < m 2 < 0.010 GeV 2 /c 4 , 
K 27T : - 0.0025 GeV 2 /c 4 < m 2 < 0.001 GeV 2 /c 4 , 

in order to distinguish between events that are consistent with a three body or a two body 
decay, see Fig. [31 

4.3 Additional selection criteria for K e3 candidates 

The additional requirements for this channel are: (1) the electron identification is 
performed by combining the calorimeter energy measurement (E) with the spectrometer 
momentum (p), and requiring E/pc > 0.95. The average electron identification efficiency 
is (98.59±0.09)%, and the efficiency as a function of p is shown in Fig. HI The error shown 
for the electron identification is only statistical in nature, but the systematic uncertainty 
is expected to be small since the result is reproduced using Dalitz decays; (2) a minimum 
momentum of 5GeV/c is required, in order to avoid low particle identification efficiency, 
and a maximum momentum of 35 GeV/c is also applied to avoid regions of very low rate; 

(3) the electron and the 7r° energies in the kaon center of mass are required to be less 
than 0.22 GeV and 0.27 GeV respectively, to further reduce the background from K 2n ; and 

(4) m^o < 0.425 GeV/c 2 is also required in order to further reduce the ir^n background. 

Using Monte Carlo simulations, the backgrounds to this channel are found to be 
well below the per mille level, see Table [2] and Fig. [31(a). The main contribution comes 
from decays with a 7r° and a tt^ that pass the E/pc requirement; for example, from K 27T 
that has a branching fraction that is about four times larger, and A' 7r ± 7r o 7r o where only the 
information from two of the photons are used. 



4.4 Additional selection criteria for candidates 

The additional requirements for this channel are: (1) in order to identify a muon 
a charged track, reconstructed in the spectrometer, is extrapolated to the MUC planes 
and associated with the MUC hits after testing the spatial separation and time difference 
between this extrapolated track and the MUC hits. Multiple scattering is taken into 
account before applying the spatial cut, and light propagation along the MUC strips is 
taken into account before applying the time-difference cut. A track is classified as a muon 
if the following conditions are satisfied: 

— each track has hits in planes one and two of the MUC; 

— for each track, the time difference measured by the chambers or the hodoscope, and 
by the muon detector does not differ by more than 4.5 ns; 

Using these requirements, the efficiency for muon identification is measured from a kine- 
matically selected sample of K — > fiu and found to be larger than 0.995 as shown in 
Fig. [51 The average muon identification efficiency is (99.759±0.003)%. The error shown 
for the muon identification is only statistical in nature, but the systematic uncertainty 
is expected to be small since the result is reproduced using data taken with a dedicated 
muon beam; (2) The individual track momenta is required to be above 10 GeV/c and be- 
low 40 GeV/c; (3) the energy of the muon in the kaon center of mass is below 0.23 GeV; 
(4) the 7T° energy in the kaon center of mass is below 0.23 GeV or m M7r o < 0.38 GeV/c 2 ; 
and (5) the 7r° energy is smaller than 40 GeV to reduce the background from K 2n . 

Backgrounds to this channel are found to be at the two per mille level, see Table [2] 
and Fig. [3](b). The main contribution comes from decays with a 7r° together with a 7r =l= 
that decays in flight. 

4.5 Additional selection criteria for K 27T candidates 

The additional requirements for this channel are: (1) muons are not rejected to avoid 
losing events where the pion decays (this corresponds to 1.6% of the K 2lT events), but 
electrons are rejected by requiring E/pc < 0.95. The average pion identification efficiency 
is (99.524±0.009)%. The fraction of surviving pions after this requirement is shown in 
Fig. H](b). The error shown for the pion identification is only statistical in nature, but 
the systematic uncertainty is expected to be small since the result at low momenta is 
reproduced using an independent sample of 7r ± 7r ± 7r =F decays. (2) The individual track 
momenta are required to be above 10 GeV/c and below 50 GeV/c. 

Backgrounds to this channel are found to be at the few per mille level, see Table |2j 
The main contribution is from because we do not reject events with muons. 

4.6 Final data samples 

Around 56K (3 IK), 49K (28K) and 462K (256K) events are found for K+ (K~ 3 ), K+ 3 
(K~ 3 ) and (K^) candidates respectively, after all the event selection requirements are 
applied (see Tabled]). 

5 Acceptance calculations 

To obtain the acceptance and an estimate of the background fractions, a GEANT[8J 
based Monte Carlo simulation is modeled for each beam separately. The beam optics for 
the K + and the K~ beamlines are described using TURTLE [9], and fine tuned using 
fully reconstructed 7r ± 7r ± 7r =F and ^tt events. 

The main inputs needed to describe the decay amplitude for semileptonic decays 
are: (1) the radiative corrections, and (2) the parameters, A+, A" and A , of the model 



presented in Eq. §E§ and used to describe the dependence of the form factors, /+ and fo, 
on the momentum transferred to the leptons, t; see Section [2J 

In order to reduce the uncertainties due to radiative effects, after identifying the 
track and the 7r°, events with extra photons are kept. Nevertheless, it is important that 
the Monte Carlo includes a well defined fraction of radiative decays (like K e ^) in order to 
obtain an acceptance correction that is accurate at the percent level. For example, studies 
show that ignoring the radiative events in the K e3 generation causes the acceptance to be 
overestimated by 1.6%, even after correcting for virtual effects on the three-body Dalitz 
plane region using the prescription given in [5]. 

The PHOTOS package [10] is used to simulate bremsstrahlung, and calculations 
from [5] are added to include virtual photons and electrons. Figures [6](a) shows a compar- 
ison between data and Monte Carlo events for the electron energy in the center of mass of 
the kaon. Good agreement is found, even though this distribution is sensitive to radiative 
effects. All other kinematic distributions show that the data are well described by the 
Monte Carlo after including effects introduced by K e3l events. For example, the invariant 
mass of the electron and ir° pair is sensitive to K e ^, see Fig. [71 The proper description 
of the low mass region proved to be sensitive to the inclusion of the K e3l . The effect 
of radiative corrections to the K e3 acceptance is confirmed using the program described 
in reference [11] that includes virtual corrections and the corresponding radiative decay 
simultaneously. The K+(K~ 3 ) acceptance is found to be 0.0709±0.0001 (0.0706±0.0001). 

For the prescription given in [5] is sufficient due to the smaller bremsstrahlung 
contribution, and adding PHOTOS had a negligible impact on the acceptance. Fig. E](b) 
shows a comparison between data and Monte Carlo events for the muon energy in the cen- 
ter of mass of the kaon. The K^ 3 (K~ 3 ) acceptance is found to be 0.0930±0.0001 (0.0927± 
0.0001). 

For a similar reason, the change in acceptance in is found to be smaller than 
one per mille if the radiative events are included in the generator. The K^iK^-n) accep- 
tance is found to be 0.1424 ± 0.0001 (0.1419 ± 0.0001). 

The response of the LKr detector, even though simulated, is not used for particle 
identification during the acceptance calculations. The particle identification efficiency is 
studied using data and corrections are made separately, see Fig. @] and [51 The acceptances 
after including the corrections due to particle identification for all three decay modes are 
summarized in Table GO, and the assigned uncertainties are statistical in nature. 

The acceptance for K + is systematically higher than for K~ due to differences 
introduced by the polarity selected for the spectrometer magnet. 

5.1 Form factor description in the Monte Carlo generation 

The Dalitz distributions for the semileptonic decays are generated with: A+ = 
0.02485 ± 0.00163 ± 0.00034, A'| = 0.00192 ± 0.00062 ± 0.00071 (both from K e3 de- 
cays) [5j, and Ao = 0.0196 ± 0.0012 (assuming /i — e universality) [2]. The change in 
acceptance when the current values of A' + , A" and A are changed by one sigma is 
about 0.08% (0.07%) for K e3 (K^ 3 ). The correlation between A + and A + is assumed to 
be -0.95 [2], while the measurements for the vector and the scalar form factor param- 
eters are assumed to be uncorrelated. Table [3] shows the the resulting slope in the 
TlKe3/K2n, T^-k^3/K2tt and TZK^/Ke?, measurements after changing form factor param- 
eter by ± one sigma, while maintaining the other parameters constant. The contribution 
from this change in acceptance to the total systematic error is included in Table HJ 

The changes in acceptance due to alternative models describing the t dependence 



of the form factors are given below: 
— the linear approximation: 



/ +>0 (i)=/ + (0) l + A +>0 — , (9) 

where A+ = 0.0296 ± 0.0008 and A = 0.0196 ± 0.0012 (assuming fi - e universal- 
ity) [2]. In this case, the ratio between the acceptance obtained using this approx- 
imation compared to the acceptance using the quadratic approximation (Eq. [5]) is 
found to be 1.0006 and 0.9998 for K e3 and K^ 3 , respectively. 
— the pole approximation: 

/+.<>(*) = /+(«) ( J% I , (io) 

\ " L v,s 

where m v = 0.877 ± 0.005 GeV/c 2 for K® 3 and K® 3 (assuming fj, — e universality) [2] 
and m s = 1.187 ± 0.050 GeV/c 2 for K® 3 (assuming fj, — e universality) [2]. These 
pole masses are obtained from neutral kaons. In this case, the ratio between the 
acceptance obtained using this approximation compared to the acceptance using 
the quadratic approximation (Eq. [5]) is found to be 0.9996 and 0.9984 for K e3 and 
K^, respectively. 

The full difference in acceptance between the pole approximation and the quadratic ap- 
proximation is assigned as an uncertainty due to the choice of form factor model, see 
Table H 

6 Result 

Table [1] lists all quantities needed to evaluate TZ^/k '■ 

Acc Kj x t tracklD , x Trig Ki x N K x (1 + A K ) 

f J^ 3 J __ (11) 

Ki/Kj Acc Ki x e tracklDi x Trig Ki x N Kj x (1 + A K .) ' 

where i,j = £3, 2n. The correction to K e3 due to particle identification efficiency (e tracklD ) 
amounts to a few percent, while the background correction (1 + A^.) is negligible. K^ 3 
and K 2n require a particle identification correction that is below the percent level, and a 
correction for background correction that is only a few per mille. 
The results for K + and K~ combined are: 

K Ke3/K27T = 0.2470 ± 0.0009(stai) ± 0.0004(syst), (12) 
T^k^Iki-k = 0.1637 ± 0.0006(staf) ± 0.0003(syst), (13) 
K Kl a/Kea = 0.663 ± 0.003(stat) ± O.OOl(syst). (14) 

The individual K + and K~ results and their systematic uncertainties are listed in Table HJ 
The sources of systematic uncertainties are due to the corrections listed in Table [1] and 
the treatment of the form factors as discussed in Section 15.11 

The ratios are found to be insensitive to the photon reconstruction and track-finding. 
Analysis of these ratios as a function of their basic distributions show stability. 

The final results are shown in Fig. [8] for lZKe3/K2ir and TZk^i,/K2-k^ and in Fig. [9] for 
Ti-K^/Kez- These can be compared to the current PDG values of TZRe3/K2n = 0.238 ± 



0.004(2], n K ^/ K2 n = 0.159 ± 0.003 [12] and K Kfl 3/Ke3 = 0.668 ± 0.008 [2J. Taking the 
current PDG value for the K 27T branching fraction, 0.2092 ± 0.0012 [2], the branching 
fractions for the semileptonic decays are found to be: 

Br(K e3 ) = 0.05168 ± 0.00019(stat) ± 0.00008(sysi) ± 0.00030(norm), (15) 
Br(K^) = 0.03425 ±0.00013(stat)±0.00006(syst)±0.00020(norm). (16) 

The uncertainty is dominated by the existing data for the K 27T branching fraction. Recall 
the corresponding PDG values [2 J are Br(K e3 ) = 0.0498 ±0.0007 and Br(K^) = 0.0332 ± 
0.0006. Higher branching fractions are found for both K e3 and K^, confirming the K e3 
results reported by the BNL-E865 collaboration [T3] . 

7 Discussion 

7.1 V us matrix element 

As discussed in Section [2j the measured partial widths for semileptonic decays can be 
combined with theoretical corrections and other experimental measurements to calculate 
\V US \. 

Using the newly measured branching fractions given in Eqs. (fT5])-(jl6]h the K± life- 
time t^ g = (1.2385 ± 0.0024) x 10" 8 s [2], Gp = (1.16637 ± 0.00001) x 10" 5 GeV" 2 [H], 
m K = 0.493677 ± 0.000016 GeV/c 2 0, S EW = 1.0230 ± 0.0003 [15J, and the full phase 
space integrals and long-distance corrections as given in Table El the |V^ S | matrix element 
times the vector form factor /+(0) is found to be (see Eq. ([3])): 

|K s |/ + (0) = 0.2193 ±0.0012, [K e3 ] (17) 

= 0.21928 ± 0.00039(stat) ± 0.00016(syst) ± 0.00062(norm) ± 0.00095(ext), 
= 0.2177 ±0.0013, [K,*] (18) 

= 0.21774 ± 0.00041(stat) ± 0.00019(s?/st) ± 0.00064(norm) ± 0.00103(ext), 

from K e3 and K^, respectively. 53% (60%) of the external (ext) error is due to the sum 
of the long-distance corrections, SU(2) and EM corrections, taken from jH El [6] for K e3 
and i^ M 3, respectively. The phase space integral, Eq. (T4]), is evaluated using the quadratic 
approximation, Eq. ([5]), and gives the next largest contribution to the external error. The 
last significant uncertainty comes from the error in the lifetime. Combining these 
| | /.(-(O) values by assuming \i — e universality, we obtain: 

|Kj/+(0) = 0.2188 ±0.0012, (19) 
\V US \ = 0.2277 ± 0.0013 (other) ±0.0019 (theo), 

where u theo" refers to the theoretical uncertainty due to /+(0), and "other" refers to 
all the uncertainties already included in Eqs. fll7l) - (|T8l) and their correlation. To extract 
| |, it is assumed that the value of /+(0) is 0.961 ± 0.008 [lj as calculated for neutral 
kaons. There is no need to make a distinction between charged and neutral kaons in /+(0) 
because the SU{2) correction is applied directly, see Eq. ©-dl]) and Table El 

These |\4s|/+(0) and \V US \ values are to be compared to predictions obtained by 
imposing the unitary condition on the CKM matrix, and taking the latest values of \ V u d\ = 
0.9738 ±0.0003 [16] and \V ub \ = (3.6 ± 0.7) x 10" 3 [2]: 

\V us \ witary f + (0) = 0.2185 ±0.0022, (20) 

\V us \ uni tary = 0.2274 ± 0.0013. (21) 



As shown in Fig. [101 these theoretical predictions are in good agreement with the results 
presented in Eq. ( TT9"j) . The result is in good agreement with the unitarity of the CKM 
mass-mixing matrix. 



The consistency between the width ratio for the semileptonic decays already pre- 
sented in Eq. (I14p and the world average for the linear form factors parameters, A + and 
A [2], can be tested. The band in Fig. [9] corresponds to the predictions for 71x^3/ Ke3 
assuming /i — e universality, Eq. ([6]), with the A + and Ao values given for K in the 
PDG of 2006 [2]. Using Eq. (ED, TZ Kl , 3/Ke3 is predicted to be 0.6682 ± 0.0017. The result 
for 'R-K V L3jKe3 suggests a lower value for Ao than the current world average for K ± [2], as 
found in recent measurements from Ki [2] decays. 

li — e universality can be tested using the TZKfj,3/Ke3 ratio and Eq. ([7]). The result 
shown in Eq. (fT4"l) implies that g^f+(0)/g e f+(0) is 0.99 ± 0.01, which is consistent with 
unity within the experimental error. 

Conclusion 

The measured ratios of the decay rates for the charged kaon semileptonic decays 
are found to be: 



Using the current experimental knowledge of the branching ratio, this leads to 
branching fractions of Br(K e3 ) = 0.05168 ± 0.00036 and Br(K^ 3 ) = 0.03425 ± 0.00024. 
This exceeds the PDG value [2] in both cases. 

Combining these results, we find |\4s|/+(0) = 0.2188 ± 0.0012, in good agreement 
with the CKM unitary prediction. 

T^Kii3/Ke3 is in reasonable agreement with the semi-empirical predictions and the 
error is improved compared to the current world average. 
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7.2 T(K, 3 )/T(K e3 ) 




0.2470 ± 0.0009(stot) ± 0.0004(s?/st), 
0.1637 ± 0.0006(stat) ± 0.0003(syst), 
0.663 ± 0.003(stat) ± 0.001(s^). 



(22) 
(23) 
(24) 



Decay type 


Raw number 
of events 


Acceptance x 
particle ID 

(AcCi X ttrack ID ) 


Backgrounds/Signal 


Trigger 
Efficiency 

(Trigi) 


Ks 


56,196 
30,898 


0.0698 ±0.0001 
0.0694 ±0.0001 


(0.0200 ±0.0008)% 
(0.0209 ±0.0010)% 


0.9990 ±0.0005 
0.9982 ±0.0008 




49,364 
27,525 


0.0927 ±0.0001 
0.0925 ±0.0001 


(0.2215 ±0.0079)% 
(0.2175 ±0.0077)% 


0.9986 ±0.0006 
0.9988 ±0.0007 


K& 


461,837 
256,619 


0.1418 ±0.0001 
0.1412 ±0.0001 


(0.2893 ±0.0058)% 
(0.2896 ±0.0058)% 


0.9987 ±0.0002 
0.9990 ±0.0002 



Table 1: Summary of information used to extract the branching ratio, where track = 
e*^, 7r ± for< = K%K%K±. 



Contributing channel 


K + 


K~ 


K e3 


K n ± n o n o 


(0.0130 ±0.0007)% 


(0.0139 ±0.0009)% 


K 2n 


(0.0070 ±0.0003)% 


(0.0071 ± 0.0004)% 






(0.1598 ±0.0045)% 


(0.1599 ±0.0046)% 




(0.0617 ±0.0064)% 


(0.0576 ±0.0061)% 


K 27r 




(0.2848 ±0.0058)% 


(0.2846 ±0.0058)% 


K e3 


(0.0045 ±0.0006)% 


(0.0050 ±0.0008)% 



Table 2: Percentage of background expected from Monte Carlo simulation for K e3 , 
and K 2n from the main contributors to their total background. 



Parameter Changed 


d'I^Ke3/K2n 


d'R'K(ja/K2ir 


d'I^Kn3/Ke3 








a; 
a; 

Ao 


-0.160 
-0.319 


-0.076 
-0.204 
-0.048 


0.113 
0.021 
-0.192 



Table 3: Dependence of the TtKe3/K2m TZk^3/K2tt an d ^-K^3/Ke3 measurements, when 
changing the form factor parameter, Aj, by ± one sigma, and maintaining the other form 
factor parameters constant. 





K 


1/K2TT 

K 




i/K2ir 

K 




3/E"e3 

K 


Central value 


0.2476 


0.2460 


0.1636 


0.1639 


0.6605 


0.6661 


Statistical error 


0.0011 


0.0015 


0.0008 


0.0010 


0.0040 


0.0055 


Total systematic error 


0.0005 


0.0006 


0.0004 


0.0004 


0.0016 


0.0018 


Accept, x Part- ID num. 
Accept, x Part-ID denom. 
Trigger efficiency in num. 
Trigger efficiency in denom. 
Background subtraction 
Uncertainty in form factor 
Uncertainty in f. f. model 


0.00041 
0.00021 
0.00013 
0.00005 
0.00001 

o.oc 
o.oc 


0.00047 
0.00021 
0.00020 
0.00005 
0.00001 

010 

009 


0.00017 
0.00014 
0.00010 
0.00003 
0.00002 
O.OC 
O.OC 


0.00017 
0.00014 
0.00011 
0.00003 
0.00002 

010 

026 


0.00067 
0.00106 
0.00039 
0.00033 
0.00005 
O.OC 
O.OC 


0.00068 
0.00126 
0.00046 
0.00053 
0.00005 

029 

086 



Table 4: Summary of the statistical and systematic uncertainties for the TZKe3/K2 
TIk^/K2it and TZx^/Kea measurements. 



Decay 


Branching 


Phase Space 


Radiative 


|K.|/+(o) 


Channel 


Fraction 


Integral 


Correctic 


njiElE] 






Br 




<W(2)(%) 


4m(%) 




K e3 


0.0517 ±0.0004 


0.1591 ±0.0012 


2.31 ±0.22 


0.03 ±0.10 


0.2204 ±0.0012 


K, 3 


0.0343 ± 0.0002 


0.1066 ±0.0008 


2.31 ±0.22 


0.20 ±0.20 


0.2177 ±0.0013 



Table 5: Inputs to Eq. (jHJ) and results for |V^ s |/+(0). By assuming unitarity, the prediction 
of |K s |/+(0) for charged kaons is 0.2185 ± 0.0022. 



n 40000 

R 

aj 

W 30000 
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2 
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10000 - 





4000 
2000 I- 





3000 
2000 
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r.m.s. 2H01 ± 0.010 (MeV/c 2 ) 




evTC A 

" J\ , 


r.m.s. 21240 ± 0.00£) (MeV/c ) 


r.m.s. 2l223 ± 0.01 1 (MeV/c ) 


- H + v/ / \ 


"■•A 


" ./ V 


- J V 




1 1 





0.12 0.13 0.14 0.12 0.13 0.14 



40000 °2 
30000 x 

OH 
4— ' 

R 

20000 % 

10000 Us 
o 

S 
4000 § 



2000 


3000 
2000 
1000 




yy Mass (GeV/c ) 



yyMass (GeV/c ) 



Figure 1: The three final states are characterized by different average photon energies 
and acceptance. This fact is reflected into the variation of the root mean square (r.m.s.) 
of the 77 invariant mass distribution. The solid curves represent Gaussian fits with the 
particular r.m.s. value noted in each plot. These variations however are small enough and 
do not introduce systematics in the normalisation of the results. A clear 7r° signal is shown 
in all channels after applying the full event selection. The factor of 1.78 corresponds to 
the ratio of K + to K~ flux in the beamline. 




0.25 0.3 0.35 0.4 0.45 0.5 0.55 

m ra0 ( GeV/ C 2 ) 

Figure 2: Monte Carlo simulation for K e %, K^, and K2 n - These are the distributions of 
the reconstructed invariant mass of the two best photons and the track without applying 
particle identification, and under the assumption that the two photons come from a 7r° 
decay and the track is a tt^. 




Missing Mass Squared (GeV 2 /c 4 ) Missing Mass Squared (GeV 2 /c 4 ) 

Figure 3: Data and Monte Carlo comparison for the reconstructed missing mass squared, 
m 2 v , for (a) K e3 and (b) K^ 3 candidates. Only the information from the charged lepton and 
the 7T° is used and extra photons are ignored, if present in the event. A good description is 
found after summing the Monte Carlo prediction for the signal and the various background 
components. Arrows indicate the allowed range in the event selection. The charged kaon 
energy is assumed to be 60 GeV. 
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Figure 4: The E/pc particle identification efficiency for (a) electrons and (b) charged pions 
from clean subsamples of K e ^ and decays respectively. These data samples have more 
stringent requirements on the kaon mass, the missing mass and the transverse momentum 
of the track. 
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Figure 5: Particle identification efficiency for muons measured using a sample obtained 
from K ± — > /i ± z/ decays. 
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Figure 6: Comparison of the (a) electron and (b) muon energy in the kaon rest frame for 
K e $ and Kfja events. Ratio of the Data over Monte Carlo is given in the lower plots. 
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Figure 7: (a) Comparison of the electron and the tt° invariant mass distribution for 
K e3 events. Ratio of the Data over Monte Carlo given in the lower plots for Monte Carlo 
generated with corrections to the Dalitz distributions for virtual effects [5 J and real photons 
due to bremsstrahlung [TU] (c) is in good agreement, while (b) only correcting the Dalitz 
distributions for virtual effects [5] and ignoring the K e3j contribution fails to described the 
low mass region. 
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Figure 8: TZKe3/K2n and 1Zk^/K2-k results compared to the corresponding PDG value [2]. 
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Figure 9: Hr^/Kbs results compared to KEK-246 results jTT], the corresponding PDG 
value of 2006 [2 J and to the predictions assuming jjl — e universality, Eq. (jSJ), with the A + 
and A values given for K ± in the PDG of 2006 f2J- 
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Figure 10: Comparison of the NA48 measurement of |\4 s |/ + (0) from K e3 and K^ 3 data, 
and the K e % BNL-E865 result [13] . The theoretical prediction shown is obtained assuming 
unitarity of the CKM matrix V u d, using the PDG values for V u d and V u b as input, and 
using the choice of /+(0) as described in the text. 
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